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Abstract 

Constitutive equations are derived for the viscoelastic behavior of particle-rein- 
forced rubbers at isothermal loading with finite strains. A filled rubber is thought 
of as a composite medium where inclusions with high and low concentrations of 
junctions between chains are distributed in a bulk material. The characteristic size 
of inhomogeneities is assumed to be small compared to the dimensions of a specimen 
and to substantially exceed the radius of gyration for macromolecules. The inclu- 
sions with high concentration of junctions are associated with regions of suppressed 
mobility of chains that surround isolated clusters of filler and its secondary network. 
The regions with low concentration of junctions arise during the mixing process due 
to the inhomogeneity in spatial distribution of a cross-linker. With reference to 
the theory of transient networks, the viscoelastic response of elastomers is ascribed 
to the thermally activated processes of breakage and reformation of strands in the 
domains with low concentration of junctions. Stress-strain relations for particle- 
reinforced rubbers are developed by using the laws of thermodynamics. Adjustable 
parameters in the constitutive equations are found by fitting experimental data in 
tensile relaxation tests for several grades of unfilled and carbon black (CB) filled 
rubber. It is evidenced that even at moderate finite deformations (with axial elon- 
gations up to 100 %), the characteristic rate of relaxation is noticeably affected by 
strains. Unlike glassy polymers, where the growth of longitudinal strain results in 
an increase in the rate of relaxation, the growth of the elongation ratio for nat- 
ural rubber (unfilled or CB reinforced) implies a decrease in the relaxation rate, 
which may be explained by partial crystallization of chains in the regions with low 
concentration of junctions. 



1 Introduction 



This paper is concerned with the viscoelastic behavior of unfilled and particle-reinforced 
rubbers at isothermal loading with finite strains. The time-dependent response of rubbery 
polymers has attracted substantial attention in the past decade, see, e.g., [^-[^, which 
may be explained by numerous applications of rubber-like materials in industry (vehicle 
tires, shock absorbers, earthquake bearings, seals, fiexible joints, solid propellants, etc.). 

Numerous observations in uniaxial and biaxial tests evidence that rubbers reinforced 
with carbon black or silica, as well as solid propellants filled with micrometer-size hard 
particles reveal strong time-dependent response which may be associated with their vis- 
coelastic behavior |10|, 0, 0, |15|, |19|, |2l[]. A mechanism for stress relaxation at 
the molecular level remains, however, an unresolved problem. 

Aksel and Hiibner |^ carried our tensile relaxation tests at ambient temperature on 
polybutadiene rubber filled with glass beads. Analyzing compounds with various con- 
centrations of filler, they reported that (i) relaxation was not observed in the unfilled 
rubber and (ii) the strength of relaxation (measured as the ratio of the increment of the 
longitudinal stress, Aa(t) = cr(0) — o-(t), at a fixed time t > to the stress, cr(0), at the 
beginning of the tests) increased with the content of reinforcement. As a mechanism for 
stress relaxation, they mentioned slippage of polymeric chains along the filler surfaces in 
the linear regime (at relatively small strains) and detachment of chains from the surfaces 
of beads and creation of vacuoles at large deformations. 



Liang et al. p2| presented experimental data in dynamic tensile tests (frequency 1 Hz) 



on ternary composite of polypropylene and EPDM (ethylene propylene diene monomers) 
elastomer filled with glass beads. They demonstrated that the mechanical damping (es- 
timated by means of the loss modulus) of the compound with untreated beads decreased 
with the filler content, in contrast with the observations of Ref. ||^. 

Leblanc and Cartault measured dynamic shear response of uncured styrene- 



butadiene rubber filled with various amounts of carbon black and silica and reported 
that the compounds exhibited rather limited time-dependent response. They concluded 
that cross-links played the key role in the viscoelastic behavior of rubbers, whereas the 
effect of filler is relatively weak. 

Clarke et al. |]T8[ carried out long-term (with the duration over 10^ s) tensile relaxation 
experiments on two unfilled elastomers (natural polyisoprene rubber with polysulphur 
cross-links and synthetic polyacrylate rubber with carbon-carbon and carbon-oxygen 
cross-links) at an elevated temperature (T = 80 °C). They reported that at relatively 
small strains (e = 0.05), natural rubber demonstrated significant mechanical relaxation, 
whereas the time-dependent response of the synthetic elastomer was not revealed. Slow 
relaxation in natural rubber was ascribed to isomerisation of polysulphur cross-links, but 
no explanation was provided for the stress relaxation in natural rubber in short-term 
mechanical tests, as well as for its absence in the polyacrylate elastomer. 

The objective of the present study is to shed some light on the mechanism of stress 
relaxation in unfilled and carbon black (CB) filled natural rubbers. For this purpose, 
we carry our tensile relaxation tests at room temperature on specimens with various 
contents of filler and develop constitutive equations for the time-dependent response of 
particle- reinforced elastomers. Adjustable parameters in the stress-strain relations are 
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found by fitting observations at various elongation ratios, A, in the region of moderate 
finite deformations, A G (1,2). 

To describe the time-dependent behavior of rubbers, we employ the concept of tempo- 
rary networks. According to the statistical theory of macromolecules, a rubbery polymer 
is thought of as a network of long chains bridged by junctions (chemical and physical 
cross-links, entanglements and van der Waals forces). A strand between two neighboring 
junctions is modelled as a set of mers connected in sequel. The classical theory of rubber 
elasticity treats all junctions as permanent, which implies that strands cannot be de- 
tached from the junctions. To predict the time-dependent response of rubbery polymers. 
Green and Tobolsky introduced a concept of temporary junctions and presumed that 
chains can slip from the junctions as they are thermally agitated. This theory was revised 
by Yamamoto Lodge and Tanaka and Edwards and was applied to the de- 
scription of the viscoelastic behavior of rubbers by Drozdov |^ and Ernst and Septanika 

According to the Green-Tobolsky concept, a strand whose ends are linked to separate 
neighboring junctions is treated as an active one. Snapping of an end of a strand from a 
junction is thought of as its breakage (transition from the active state of the strand to its 
dangling state). When a dangling chain captures a junction (which reflects creation of a 
new physical cross-link or entanglement), a new active strand merges with the network. 
Breakage and reformation of active strands are treated as thermally activated processes: 
attachment and detachment events occur at random times as they are driven by thermal 
fluctuations. 

A particle-reinforced rubber is modelled as a micro-composite, where domains with low 
and high concentrations of junctions are distributed in a matrix with a "normal" concen- 
tration of cross-links between chains. Regions with low concentration of junctions (RLCJ) 
arise during the mixing process because of the inhomogeneity in the spatial distribution 
of a cross-linker. Regions with high concentration of junctions (RHCJ) are formed in a 
particle-reinforced elastomer in neighborhoods of filler clusters, where the filler particles 
and their aggregates play the role of extra cross-links between polymeric chains. Breakage 
and reformation of strands are assumed to take place in the regions with low concentra- 
tion of junctions only, whereas in the elastomeric matrix and in the domains with high 
concentration of junctions slippage of chains from cross-links is prevented by surrounding 
macromolecules. The viscoelastic response is ascribed to the reformation of strands in 
RLCJs only, whereas the host matrix is treated as a nonlinear elastic medium. We pos- 
tulate that the response of domains with high concentration of junctions is substantially 
affected by mechanical loading, which results in damage (rupture) of isolated clusters of 
filler and its secondary network. Mechanically induced changes in the topology of filler 
clusters are described by a vector of internal variables, a, which is not concretized in the 
present study. To distinguish between the time-dependent response of the composite in 
the standard relaxation tests driven by reformation of strands in RLCJs and that induced 
by changes in the internal variables, we presume that the derivative of the vector a with 
respect to time vanishes when the rate-of-strain tensor equals zero, see Eq. (48). 

The aim of this paper is to derive a constitutive model for the time-dependent behavior 
of particle-reinforced rubbers and to apply this model to the analysis of experimental 
data in tensile relaxation tests. Our purpose is to assess the effects of filler content and 
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elongation ratio on the relaxation spectrum and the rate of stress relaxation. 

The study is organized as follows. The free energy density of a particle-reinforced 
elastomer is determined in Section 2. Stress-strain relations at finite strains are derived 
in Section 3 by using the laws of thermodynamics. Uniaxial extension of a specimen is 
analyzed in Section 4. Section 5 deals with experimental data in tensile relaxation tests 
for three different grades of CB filled rubber. Adjustable parameters in the constitutive 
equations are found in Section 6 by fitting observations. A brief discussion of our findings 
is presented in Section 7. Some concluding remarks are formulated in Section 8. 



2 Free energy of particle-reinforced rubber 

A particle-reinforced rubber is modelled as a micro-composite where inclusions with high 
and low concentrations of junctions are randomly distributed in the bulk material. The 
composite is treated as an incompressible and isotropic elastic medium with the free 
energy density (per unit mass) 

XiJ = + V + W, (1) 

where \E'° is the free energy in the stress-free state at the reference temperature T°, V is 
the free energy of thermal motion of chains, and W is the increment of free energy driven 
by mechanical deformations. 

We confine ourselves to quasi-static loadings with moderate finite strains, when me- 



chanically induced changes in temperature, AT = T — T°, are rather weak |31|. This 
implies that the effect of temperature on the specific heat, c, as well as thermal expansion 
of the elastomer may be neglected. Assuming c to be temperature-independent, we adopt 
the conventional formula for the free energy of thermal motion 

V^=(c-5°)(T-n-cTln^, (2) 

where 5"° is the entropy per unit mass in the stress-free state at the reference temperature. 

Let /iL and /xh be mass fractions of regions with low and high concentrations of junc- 
tions. Taking into account that the strain energy density, W, is an extensive thermody- 
namic potential, we find that 

W={l-flL-m)WB+f^LWL + mWH, (3) 

where Wb is the mechanical energy per unit mass of the bulk medium, and PFl, are 
the strain energy densities per unit mass for regions with low and high concentrations of 
junctions, respectively. 



2.1 Mechanical energy of the host matrix 

The bulk medium is treated as a permanent network of long chains bridged by junc- 
tions (chemical cross-links, physical cross-links whose life-time exceeds the duration of a 
mechanical test, and entanglements). We adopt the affinity hypothesis which neglects 
thermal oscillations of junctions and postulates that the displacements of junctions at 
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the micro-scale coincide with the displacements of appropriate points of a sample at the 
macro-level. 

To clarify this assumption, we consider a strand between two contiguous junctions 
Ai and In the stress-free state, the junctions occupy points with the radius vectors 
R° and Denote by Ri(t) and R2(t) radius vectors of these junctions at an arbitrary 
instant t > (the time t = corresponds to the instant when external loads are applied to 
the specimen). Let r° = Rg — R^ be the end-to-end vector for the strand in the stress- free 
state and r(t) = R2(t) — Ri(t) the end-to-end vector in the deformed state at time t. 
According to the affinity hypothesis, the deformation gradient Vr°r(t) for motion of junc- 
tions coincides with the deformation gradient VRoR(t) = V°R(t) for macro-deformation. 
This implies that the Cauchy deformation tensor c{t) for the network equals the Cauchy 
deformation tensor for macro-deformations, 

C°(t) = V°R(t) ■ [v°R(t) 

and the rate-of-strain tensor for the network d(t) coincides with the rate-of-strain tensor 
for macro-deformations 



D(t) = ^[V(t)V(t) + (v(t)V(t))^]. (4) 
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Here V(t) = VR(t) is the gradient operator in the deformed state of the specimen, 

is the velocity vector, the dot stands for inner product, and T denotes transpose. 

For an isotropic incompressible medium, the strain energy density Wb depends on the 
current temperature T (the strong effect of temperature reflects the entropic nature of 
the mechanical energy for elastomers) and the flrst two principal invariants of the Cauchy 
deformation tensor C°, 

W^it) = WB{Tit),hiC%t)),hiC°it))). (5) 

Disregarding thermally induced changes in the reference state (at small increments of 
temperature, AT, of the order of several K ||3l]] thermal expansion of rubber is negligible 
compared to the strains under consideration), we assume the function to obey the 
condition 

WB{T,h,h)l^^^^^^^=0, (6) 

which means that the strain energy vanishes in the stress-free state. 
The mechanical energy per strand of the bulk material is given by 

WB{t) = (7) 
where A^b is the average number of strands per unit mass. 
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2.2 Mechanical energy of domains with low concentration of 
junctions 

Regions with low concentration of junctions are assumed to be randomly distributed in 
the host matrix as inclusions with various sizes and shapes. Any RLCJ is modelled as a 
transient network of long chains, where breakage and reformation events occur at random 
times as active and dangling strands are agitated by thermal fluctuations. With reference 
to the theory of thermally activated processes , we assume that the rate of breakage 
for a strand in a stress-free medium is given by 



roexp(- 



where a) is the energy for breakage, is Boltzmann's constant, and the pre-factor Fq 
weakly depends on temperature. It follows from Eq. (8) that for any > 0, F vanishes at 
low temperatures, T — 0, and F tends to Fq at high temperatures, T — > oo. This implies 
that the coefficient Fq may be associated with the breakage rate at elevated temperatures. 
We introduce the dimensionless energy for breakage. 



kBT°' 

and neglect the effect of small increments of temperature, AT, on the rate of breakage. 
This results in the formula for the rate of transition from the active state of a strand to 
its dangling state 

F(cu) = Foexp(-c^). (9) 

Equation (9) is assumed to be satisfied for an arbitrary loading, provided that the coeffi- 
cient F depends on the Cauchy deformation measure, C°, and additional internal variables 
that characterize damage of the secondary network and isolated clusters of filler. Because 
a detailed description of damage for particle-reinforced composites is beyond the scope of 
the present study, we accept a formal standpoint and postulate that rupture of filler clus- 
ters is described by a vector of internal parameters a. This implies that for an isotropic 
and incompressible medium, 

Fo = Fo(a,/i(C°),/2(C°)), (10) 

where the function Fo(a, Ji,J2) is to be found by fitting experimental data. It follows 
from Eqs. (9) and (10) that 

F(t,w) = F,(t)exp(-cu) (11) 

with 

F,(t)=Fo(a(t),Ji(C°(t)),J2(C°(t))). 

Denote by p{uj) the distribution function for energies of breakage for strands in RL- 
CJs. The function p{u) is assumed to be robust with respect to the infiuence of thermo- 
mechanical factors. This means that p is treated as a function which is temperature- 
independent (for small increments of temperature, AT) and which does not depend on 
the history of loading. To fit experimental data, we use the quasi-Gaussian distribution 



p{uj) = po exp 



2n 



2S2 
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where fl and S are adjustable parameters, and the constant po is found from the condition 

/ p(uj)du} = 1. (13) 
Jo 

Let dN{uj) be the concentration of strands in RLCJs (per unit mass) whose energy for 
breakage, a;', is located in the interval [a;, a; + duj]. The quantity dNiuj) reads 

dN{uj) = Ni,p{uj)duj, 

where A^l is the average number of strands per unit mass of RLCJs. 

The kinetics of reformation is described by the function X{t,T,uj), which equals the 
current (at time t) number (per unit mass) of active (merged to the network) strands with 
the energy for breakage uo that have last been bridged to the network before an instant 
r G [0,t]. The function X entirely determines the current state of active strands in the 
regions with low concentration of junctions. The quantity X{t,t,uj) equals the number 
of active strands (per unit mass) with the energy for breakage uj at time t. In particular, 
X(0, 0,0;) is the initial number of active strands with the energy for breakage uj, 

X(0,0,a;) = iVLP(^). (14) 

The amount 



dX. 



dr 

t=T 



is the number (per unit mass) of active strands with the energy for breakage cu that have 
been linked to the network within the interval [r, r + dr] . The quantity 

— {t,T,u) dr 

determines the number of these strands that have not been broken during the interval 
[r, t] . The amount 

- — (t,0,u;)dt 

is the number of active strands (per unit mass) that detach from the network (for the 
first time) within the interval [t,t + dt], and the quantity 

d'^X 

^ t7t~(^''^'^) ^^^^ 

OtOT 

is the number of strands (per unit mass) that have last been linked to the network within 
the interval [r, r + dr] and leave the network (for the first time after merging) during the 
interval \t,t-\- dt\. 

The kinetics of evolution for a transient network is determined by the relative rate 
of breakage for active strands, T{t,uj), and by the rate of merging with the network 
for dangling strands, 7(t, a;). The assumption that F and 7 are functions of energy for 
breakage, a;, distinguishes the present model from previous ones. The dependence of the 
rates of breakage and reformation on time, t, refiects the effects of current strains and the 
damage variables, a, on these parameters under time-dependent loading. 
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The relative rate of breakage, F, equals the ratio of the number of active strands 
broken per unit time to the total number of active strands. Applying this definition to 
the strands arisen before loading and remaining active at the current time t, we find that 

1 dX 

Applying this definition to the strands merging with the network at time r and remaining 
active at time t > r, we obtain 



dtdr 



{t,T,Uj) 



(16) 

the rate of refor- 



In agreement with conventional theories of temporary networks |^ 
mation, 7, is defined as the number of dangling strands (per unit mass) bridged to the 
network per unit time 

. N dX^ 



dr 



(17) 



We confine ourselves to a model with two possible states (active and dangling) for a strand. 
A dangling chain is treated as a chain where stresses totally relax after its detachment (no 
memory about previous deformations), whereas an active one is a chain that preserves 
the entire memory about its state at the instant when it merges with the network. One 
can hypothesize about an intermediate state between these two extremities, when the 
interval between the detachment event and the subsequent attachment to the network is 
too small for total relaxation of stresses, but is sufficiently large to ensure that a part of 
these stresses do relax [^. This is equivalent to the assumption that a partial memory 
preserves in a strand about the history of its deformation. To simplify the analysis, this 
case is excluded from our consideration. 

Formulas (15) and (16) can be treated as ordinary differential equations for the func- 
tion X. Integration of these equations with the initial conditions (14) and (17) implies 
that 



X(t,0,w) 
dX 



dr 



{t,T,Uj) 



Ni,p{uj) exp 
7(r, u) exp 



— / T(s,uj)ds 
Jo 

— J T{s,uj)ds 



;i8) 



To exclude the rate of reformation, 7, from Eq. (18), we postulate that the concentrations 
of active and dangling strands are independent of macro-strains. This means that the 
number of strands broken per unit time coincides with the number of strands merging 
with the network within the same interval. The number of strands detached from the 
network during the interval [t, t + dt] equals the sum of the number of initial strands (not 
broken until time t) that slip from temporary junctions 

- — {t,0,uj)dt 

and the number of strands linked with the network within the interval [r, r + dr] (for the 
last time before instant t) for various instants r and broken within the interval [t, t + dt] 

d^X 
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This number coincides with the number of strands, 7(^,0;) dt, merged with the network 
within the interval [t,t + dt], which results in the balance law 



'~dt 



{t,0,u;) 



dtdr 



{t, T, Uj)dT. 



Substitution of Eq. (18) into this equality implies that 



7(t,a;) = r(t,a;)<^ iVLp(a;) exp 



r(s, uj)ds 



To solve this equation, we introduce the notation 



7(r, a;) exp 



7(^,0;) = 7(t,a;)exp 



/ Vis,ijj)ds 
Jo 



and find that 



Jo 



^{t,uj) = T{t,uj) 

Setting t = in Eq. (20), we obtain 

7(0,u;) = A^Lr(0,a;)p(u;). 

Differentiation of Eq. (20) with respect to time results in 

1 dT 



r{t,u) dt 

Integration of this equality from zero to t yields 



In 



7(^,0;) 



In 



T{t,u) 



7(0,0;) r(0,a;) 
It follows from Eqs. (21) and (22) that 



+ / T(s,ijj)ds. 
Jo 



fvis,u)dt 
Jo 



r(s, u})ds 



dry 



(19) 
(20) 

(21) 



(22) 



7(^,0;) = A^Lp(<^)r(i,a;)exp 

Combining this equality with Eq. (19), we find the rate of reformation 

^(t,oj) ^ N^p(u;)r{t,u;). 



(23) 



Denote be wl the mechanical energy of an active strand in domains with low concentration 
of junctions. Because RLCJs are treated as isotropic and incompressible media, wl is 
assumed to depend on temperature, T, and the first two principal invariants of the Cauchy 
deformation tensor. The conventional hypothesis that stress in a dangling strand totally 
relaxes before this strand captures a new junction implies that the reference (stress-free) 
state of a strand merging with the network at time r coincides with the deformed state of 
the network at that instant. This means that for a strand not broken within the interval 

[o,t], 

w^(t,0) ^w^{T(t),h(C°(t)),h(C°(t))), (24) 
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and for an active strand which has last been reformed at time r e [0, t] , 



w^t, t) = w^{T{t), 7i(C(i, r)), h{C{t, r))) , 



(25) 



where C(t, r) is the relative Cauchy deformation tensor for transition from the deformed 
state of a specimen at instant r (when the strand has last been merged to the network) 
to its deformed state at time t > r. The tensor C(i, r) is expressed in terms of the 
deformation gradient V(r)R(i) by the formula 

C{t,T) = V(T)R(t) • [V(T)R(t) " 

Summing mechanical energies for strands with various energies for breakage, uj, merged 
with the network at various instants, r G [0,t], and neglecting the energy of interaction 
between strands, we arrive at the strain energy density per unit mass of RLCJs 



dcu. 



(26) 



(27) 



poo r rt QX 

Wi^{t)= / X(t,0,u;)«;L(t,0)+ / —{t,T,uj)wi^(t,T)dT 

Jo I Jo UT 

Substitution of expressions (18), (23), (24) and (25) into Eq. (26) imphes that 

P OO f ft 

W^{t) = p{u)du\w^{T{t),h{C°{t)),h{C°{t)))exv[- jj{s,uj)di 
+ r(r, u)W^[T{t), h{C{t, r)), h{C{t, r 
X exp — J r{s,uj)ds o?t|, 
where, in accord with Eq. (7), we use the notation 



2.3 Mechanical energy of domains with high concentration of 
junctions 

Because the present study focuses on the viscoelastic response of filled rubbers, we do not 
dwell on changes in the mechanical energy of regions with high concentration of junctions 
driven by rupture of the secondary network and isolated clusters of filler. It is assumed 
that the strain energy density per unit mass of RHCJs reads 



Wn{t) = Wn{3L{t),T{t), /i(C°(i)), /2(C°(i))), 
where no concrete form is provided for the function iyH(ci, T, Ji, I2). 

2.4 Mechanical energy of the compound 

Substitution of expressions (5), (27) and (28) into Eq. (3) results in the formula 

W(t) = (l-/XL-//H)W-B(n^),^l(C°(t)),/2(C°(i))) 
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(28) 



rp{Lj)du;lWj^(T{t),h{C°{t)),l2{C°{t)))exp[- j T{s,uj)ds 



+ r(r,cu)W^L(T(t),/i(C(t,r)),/2(C(t,r))jexp[- / V{s,uj)ds 
+/XHW^H(a(i), T{t), h{C°{t)), l2{C\t))). 
Bearing in mind Eq. (6), we find from Eq. (29) that 



dr 



(29) 



= J(t)^(t)+A(t).^(i)+C/°(t)^(C°(t))+C/°(t)^(C°(t))+C/(t)-y(t). (30) 
Here we introduce tlie notation 



J{t) 



(1 - /^L - /XH)^(T(t), /i(C°(t)), HC^t))) 

/^~M^)dc.|^(T(i),7i(C°(i)),72(C°(i)))exp[- j\{s,uj)ds 
+ /Jr(T,u;)^(r(i),/i(C(i,T)),/2(C(i,T)))exp[- jy{s,u;)d 



dr 



dWn 



dT 



A{t) 

m) 



[a(t),T(t),/i(C°(t)),/2(C°(t))j, 

//H Val^H (a(i) , r(i) , /i (C° (i) ) , /2(C° (i))) , 

(l-/^L-/XH)^(T(t),/l(C°(t)),/2(C° 

/o (94 

+/.H^(a(t),r(t),/i(C°(t)),72(C°(t))), 



+I^L ^ (T(t), /i(C°(i)), /2(C°(t))) exp [- T{s, u)ds\p{u)du 



U{t) = /.L /^"p(a;)cia; /J r(r, a;) [^(r(i), /i(C(t, r)), /2(C(t, r))) ^(C(i, r)) 





exp 




ci;)ds 


y(i) = //L/^"r(i,a;)p(u;){iyL(T(i),7i(C°(t)),72(C°(i))) 


exp 




uj)ds 



1^ r(r, u;)I^L(r(i), /i(C(i, r)), /2(C(i, r 



+ 
X exp 



— / T{s,uj)ds 



dr (-dcu. 



(31) 



The principal invariants of the Cauchy deformation tensor, C(t,r), coincide with the 
principal invariants of the Finger tensor for transition from the deformed state at time r 
to the deformed state at time t, 



F(t,r) = V(r)R(t) ■ V(r)R(t). 
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This assertion implies that 



f^(C(t,r)) = ^(F(t,r)) = ||(F(t,r)): 



(32) 



where the colon stands for convolution of tensors. The derivatives of the principal invari- 



ants of the Finger tensor read 



dF 



I, 



OF 



(F) = /i(F)I-F' 



where I is the unit tensor. The derivative of the tensor F with respect to time is given 

by a 

^(t,r)= [v(t)V(t)]^F(t,r) + F(t,r)-V(t)V(t). (33) 
It follows from Eqs. (32) and (33) that 
dh 



^(C(t,r)) = 2F(t,r):D(t), 



dt 



;C(t,r)) 



/i(C(t,r))F(t,r)-F2(t,r) : D(t), 



(34) 



where the rate-of-strain tensor D is determined by Eq. (4). By analogy with Eq. (34), 
we find that 



dh^ 
'dt' 



C°(t)) = 2F°(t):D(t), 



^(C°(t)) = 2 



/i(C°(t))F°(t)-(F°(t) :D(t) 



(35) 



where 



F°(t) = [V°R(t)J ■V°R(t). 
Substitution of expression (34) into formula (31) for the functional U{t) implies that 



Here 



U{t) = 2/iLAL(t) : D(t). 



p{uj)dto I 0Li(t,T)F(t,r) +0L2(t,r)F'(t,r) 

rt 

dr. 



(36) 



xr(r, ci;)exp — / T{s,uj)ds 
and the functions 4>Lk{t,T) read 



(37) 



0Ll(t,T) 



dh 



{T{t),h{C{t,r)),h{C{t,r)) 



+h{C{t, T))—^{T{t), h{C{t, r)), h{C{t, T 



0L2(t,T) 



dWi^ 
' dh 



dh 

T{t),h{C{t,T)),h{C{t,T))). 
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It follows from Eqs. (31) and (35) that 

dh 
dt 



2 Jj 

T.Um-7^{C%t)) = 2{(l-/.L-/iH)ABW+/^HAHW 



k=l 



+/iLA^(t) r 

Jo 



exp 



r{s,uj)ds\p{uj)duj} : D(t), (38) 



where 



A;(t) = 0^i(t)F°(t) + 0;2W[F°(t) 
and the functions (t>°^k(t) are given by 

OWb 



(/? = B,L,H), 



0LlW 



5/l 



T(t),/i(C°(t)),/2(C°(t)) 



^i(C°W)^(T(t),/i(C°(t)),/2(C° 



5Wl 
5/i 



r(t),/i(c°(t)),/2(c°(t)) 



0HlW 



+/i(C°(t))^(T(t),/i(C°(t)),/2(C 
''''^^a(t),T(t),/i(C°(t)),/2(C°(t)) 



/i(C°(t))^(a(t),T(t),Ji(C°(t)),/2(C° 



0H2W 



(T(t),/i(C°(t)),/2(a 

T(t),/i(c°(t)),/2(c°(t)); 

(a(t),T(t),/i(C°(t)),/2(C°(t))). 



Combining Eqs. (30), (36) and (38), we arrive at the formula 



with 



^(t) = m^it) + A(t) ■ ^(t) + 2H(t) : D(t) - Yit) 



H(t) = (l-/iL-/^H)A°(t)+/iHA^(t) 

+Hi^^A'l{t) J exp - J T{s,uj)ds p{uj)duj + Ai^ 



3 Constitutive equations 



For an incompressible medium, the Clausius-Duhem inequality reads 
.dQ 



T- 



dt 



dT d^i _ 1 



(39) 



(40) 



(41) 



(42) 
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where p is mass density, S' is the deviatoric component of the Cauchy stress tensor S, q 
is the heat flux vector, S is the entropy per unit mass, and Q is the entropy production. 
It follows from Eqs. (1), (2) and (40) that 

dt ^ T° ^ dt dt ^ ^ 

Equations (42) and (43) imply that 



P 



+A-^ + y--l;q-Vr>0. (44) 

Because inequality (44) is to be satisfied for an arbitrary loading program, the expressions 
in brackets should vanish. This results in the formula for the entropy per unit mass 

S{t)^S° + c\n^-J{t) (45) 

and the constitutive equation 

S(i) = -P(i)I + 2p{(l-/XL-/XH)A°(t)+/XHA^(i) 

[Al{t) exp(- T{s, uj)ds)p{uj)dw + AL(t)] }, (46) 

where P is pressure. To ensure that the rate of entropy production is non-negative, it is 
natural to suppose that the derivative of the vector of internal variables a with respect to 
time is proportional to the vector A, which implies the evolution equation 

|=-,A. (47) 

where 77 > is a function to be determined. We assume that t] depends on the vector a, 
the rate-of-strain tensor D, and the Cauchy deformation tensor C°, 

77 = 77(a,D,C°). 

To distinguish explicitly the viscoelastic behavior of filled rubber and damage of filler 
clusters, it is postulated that the coefficient r] vanishes for a time-independent loading 
process 

77(a,D,C°)|^^^=0. (48) 

By analogy with Eq. (47) , we adopt the Fourier law for the heat fiux vector q, 

q = -AcVT (49) 

where k > is a coefficient of thermal difi^usivity. 

Substituting expressions (45) to (47) and (49) into Eq. (44) and bearing in mind 
formula (31) for the functional Y , we find that for any loading process, the Clausius- 
Duhem inequahty (44) is satisfied. 
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4 Uniaxial extension of a specimen 



We now apply constitutive equation (46) to determine stresses in a bar at uniaxial exten- 
sion. Points of the bar refer to Cartesian coordinates {Xi} in the reference (stress- free) 
state and to Cartesian coordinates {xi} in the deformed state, {i = 1, 2, 3). Tension of an 
incompressible medium is described by the formulas 



xi^k{t)Xi, X2^k 2(^)^2, xs^k 2(^)^3, 



(50) 



where k — k{t) is the extension ratio. It follows from Eq. (50) that the deformation 
gradient, V°R(t), and the relative deformation gradient, V(T)R(i), read 



V°R(i) = k{t)eiei + 



kit) 



(6262 + 6363), 



V(r)R(t) = ^ei6i+f^) '(6262 + 6363), 



\k{t)J 



where ej are base vectors of the frame {X^}. The Finger tensor, F°(i), and the relative 
Finger tensor, F(t, r), are given by 



F°{t) = A:'(t)ei6i + -^(6262 + 6363), 
F(^,t) = (^^j eiei + ^(6262 + 6363). 



k{t) 



(51) 



Substituting Eqs. (37), (39) and (51) into Eq. (46), we find the non-zero components 



P{t) + 2p\ (1 - /XL - /^h) 0Bl(^) + <pB2it)k'it) 



cxp 



T{s,uj)ds p{uj)duj \k^{t) 



E2(i) 



+2pMLf pM<^{[0L:(*,r)(M)V^..(t.r)(M)' 

xr(r, a;)exp — j T(s,uj)ds dr, 
-P{t) + 2p|(l - /XL - /xh) [0biW + M)k-\t) 

r 1 

+/^L <Pli{t) + (Pl2{t)k-\t) / exp 

L J Jo 

k{T 

Mt) 



+ m\(l>ni(t) + (l>H2(t)k-\t) 



— / r(s,uj)ds p(Lj)dLj>k ^{t) 



+2p,, fpi.yi. [ [Mt, r) (^) + Mt, r) i^^f 



xr(T, a;)exp — j T[s,ijj)ds dr 
of the Cauchy stress tensor 

S = E16161 + ^2(6262 + 6363). 



(52) 
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Excluding the unknown pressure P{t) from Eq. (52) and the boundary condition 
we arrive at the formula for the longitudinal stress 



Ei(i) = 2p|(1-/xl-/xh)0biW+/^h0hiW 



+ 



r(s, u)ds p{u})du} 



(1 - //L - )Uh)(/>B2W + /^H(/'h2 W 



X 



A;2(i) + A;-i(i)l UA;2(i) -A;-^(i) 



p{uj)duj <^0Ll(t,T) +0L2(^,t) 



X 



u(t); k{t) 



\k{r)J ^ k{t) 



r(r, a;)exp — / T{s,uj)ds 



dr. 



(53) 



Equation (53) describes the mechanical response of a specimen for an arbitrary loading 
program. To compare results of numerical simulation with experimental data, we focus 
on the relaxation test with 

' 1, t<0, 
A, t>0, 



(54) 



where A > 1 is a constant. Substitution of Eq. (54) into Eq. (53) results in 

Si(i) = [(7b + ctl + an - a{t)] (A^ - A'^), 

where the parameters 

CTb = 2p{l - fXi^- Hu) (l)°Bl+ 4>°B2i^'^ + ^~^) : 



(55) 



<7l 



2p/iL 



LI +0L2 



(A^ + A-^) 



(7h = 2p//h[0hi + </'h2(A' + A-^) 
are independent of time and 

a{t) — ai, J 1 — expf— y r{s,u!)ds 

It follows from Eqs. (13), (55) and (56) that the dimensionless ratio 

Si(t) 



p{uj)djjj. 



(56) 



Ei(0) 



is given by 



Rit) 

Jo 



m 

1 — exp^— exp(— cc;)t 
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p{uj)du!, 



(57) 



where 

A = . (58) 

0"B + O'L + CTH 

Equations (12) and (57) are determined by 4 adjustable parameters: the constant A is the 
ratio of the relaxing stress to the total stress at the beginning of the test, F^, characterizes 
the average rate of relaxation, VL and S in Eq. (12) describe the distribution of relaxation 
times. 

Our objective now is to find the constants A, F^, Q and S by fitting experimental data 
in relaxation tests with various elongation ratios. A, for unfilled and CB filled rubbers. 



5 Experimental 

To assess the time-dependent response of particle- reinforced elastomers, a series of uniaxial 
relaxation tests were carried out at a constant temperature. Dumbbell specimens were 
provided by TARRC laboratories (Hertford, UK) and were used as received. 

The first compound on the base of natural rubber (EDS-19) contains 1 phr (parts 
per hundred parts of rubber) of carbon black for UV protection. The filler content is too 
low to affect the mechanical response significantly, and we refer to this compound as an 
unfilled rubber. 

The other compound (EDS-16) contains 45 phr of carbon black and this compound 
is treated as a filled rubber. Before the vulcanization process, a sheet of rubber from 
which the specimens were cut out was milled in a rolling mill to destroy aggregates of 
filler particles. This procedure leads to a weak anisotropy of the material. The dumbbell 
specimens cut in the direction of milling are referred to as Rl, whereas the samples cut 
in the orthogonal direction are abbreviated as R2. 

A detailed chemical formulation of compounds EDS-16 and EDS-19 is presented in 
Table 1. 

The relaxation tests were performed at room temperature by using a testing machine 
designed at the Institute of Physics (Vienna, Austria). A specimen was loaded with the 
strain rate 5.0 • 10~^ (s~^) up to a given elongation ratio A, which was preserved constant 
during the relaxation test (1 hour). To measure the longitudinal strain, two reflection 
lines were drawn in the central part of each specimen before loading (with the distance 
7 mm between them). Changes in the distance between these lines were controlled by a 
video-extensometer (which ensured the accuracy of about 1 %). 

The tensile force was measured by using a standard loading cell. The nominal stress 
was determined as the ratio of the axial force to the cross-sectional area of a specimen (1 
mm X 4 mm) in the stress-free state. 

The relaxation process was studied at three elongation ratios (A = 1.2, A = 1.4 and 
A = 1.8) for samples made of the unfilled rubber, and at four elongation ratios (A = 1.2, 
A = 1.4, A = 1.8 and A = 2.0) for specimens made of the filled rubbers Rl and R2. 
None of the specimens was subjected to heat treatment or mechanical pre-loading prior 
to testing. 
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6 Comparison with experimental data 



We begin with matching observations in relaxation tests for unfilled rubber. First, we fit 
experimental data in the test with the minimum elongation ratio Amin = 1-2. Because 
the parameters F^, and Q appear to be mutually dependent (an increase in Q leads to 
the growth of F^,), we set F* = 1.0 (s~^) in the approximation of experimental data at 
Amin- Given A, the quantities Q and S are found by the steepest-descent method. The 
coefficient A is determined by the least-squares technique. 

Afterwards, we match observations for the unfilled rubber at the elongation ratios 
A = 1.4 and A = 1.8 by using the parameters Q and S which were found in fitting 
data at Amin- The relaxation rate F^, is determined by the steepest-descent procedure, 
and the coefficient A is found by the least-squares method. Figure 1 demonstrates fair 
agreement between observations and results of numerical simulation. It is worth noting 
that relaxation curves at different elongation ratios. A, noticeably differ from each other, 
despite some scatter in the experimental data. 

The average relaxation time, r^, = F~^, and the dimensionless coefficient A are plotted 
in Figures 2 and 3 versus the first invariant of the Cauchy deformation tensor, /i(C°). The 
effect of straining on and A is adequately described by the phenomenological equations 

T-* = ro + ri(/i-3), A = Ao + Ai{h-3), (59) 

where the constants Tk (tq > 0) and A/, {k = 0, 1) are determined by the least-squares 
method. 

We proceed with matching experimental data for two grades (Rl and R2) of filled 
rubber by applying the same procedure of fitting. Observations and results of numerical 
analysis are depicted in Figures 4 and 5 which demonstrate rather small scatter of exper- 
imental data. The parameters r* and A are plotted in Figures 6 and 7 together with their 
approximations by Eq. (59). 

To compare our findings with other observations, we apply the same procedure of 
fitting to data for a filled rubber obtained in tensile relaxation tests by Haupt and Sedlan 
pT| (the samples used in their study are abbreviated as R3). For a description of speci- 
mens and the experimental procedure, the reader is referred to the paper Figure 8 
demonstrates good agreement between observations and results of numerical analysis. The 
experimental constants and A are plotted in Figures 9 and 10 together with predictions 
ofEq. (59). 

Adjustable parameters for various grades of unfilled and CB filled rubbers are listed 
in Table 2. 



7 Discussion 

Experimental data depicted in Figures 1, 4, and 5 reveal that the response of the rubbers 
under consideration (both carbon black filled and unfilled) in the standard relaxation tests 
strongly depends on time. The relative decrease in stresses after relaxation for to = 1 hour 
at the minimum elongation ratio Amin = 1.2, 

_ Si(0)-Si(to) 
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equals 6 % for the unfilled rubber and varies from 13 to 15 % for the CB filled rubbers. The 
quantity C(^o) for the filled rubbers (Rl and R2) noticeably exceeds (about by twice) that 
for the unfilled rubber. The fact that the presence of filler accelerates relaxation in rubbers 
may be explained by an increase in the number of regions with low concentration of cross- 
links during the mixing process in a particle-reinforced elastomer compared to that for an 
unfilled medium. However, because the amounts C(^o) have the same order of magnitude 
for unfilled and filled rubbers, we suppose that the viscoelastic behavior of elastomers 
cannot be described exclusively by means of the mechanically induced destruction of CB 
clusters and de- wetting of macromolecules from the filler particles [||. 

The relaxation spectrum (which is characterized by the distribution function, p{uj), for 
strands with various energies for breakage, a;) is independent of the elongation ratio. A, at 
least, in the region of moderate deformations, A G (1, 2). The distribution function, p(a;), 
is determined by two adjustable parameters, Q and S, which accept similar values for the 
three grades of filled rubber (Rl, R2, and R3). Their values, however, substantially differ 
from those for the unfilled rubber. Because the ratio ^ = lies in the interval between 
0.5 and 0.6 for all specimens, we may associate Q with the average energy and S with the 
standard deviation of energy for breakage. Table 2 demonstrates that the average energy 
for breakage of strands, Q, for the unfilled rubber exceeds by twice (approximately) the 
average energy for breakage for the filled rubbers, whereas the ratio C, (that characterizes 
the width of the Gaussian distribution) is practically independent of the presence of filler. 
This observation means that reinforcement of rubber by CB particles results in an increase 
in the level of inhomogeneity of the compound and in the growth of the number of regions 
with low concentration of junctions where reformation of strands occurs. 

Figures 2, 6, and 9 evidence that the characteristic time of relaxation, t^,, increases 
with the growth of longitudinal elongation A. This conclusion contradicts the strain- 
time superposition principle conventionally employed in the analysis of the viscoelastic 



response of polymeric glasses, see, e.g., [Bq, 37]. According to the free-volume concept 



mechanical loading results in an increase in the free volume, which, in turn, causes an ac- 
celeration of the reformation process for mobile units (cooperatively rearranging regions). 
The free-volume theory correctly describes the time-dependent behavior of glassy poly- 
mers, where mechanically induced changes in the specific volume are substantial. The 
applicability of that approach to rubbery polymers appears to be questionable, because 
of the incompressibility of rubbers in the region of moderate strains (up to A = 3 for 
uniaxial tension). 

The growth in the characteristic time for rearrangement of strands, r^,, with A may 
be explained by the strain-induced crystallization of natural rubber. Uniaxial extension 
of a specimen at the macro-level causes stretching of strands at the micro-level. This 
results in the creation of nuclei of crystallites (regions of strong order and high density) 
in the vicinity of junctions, which slow down breakage and reformation of chains. With 
reference to this picture, one can hypothesize that the presence of filler particles (which 
may be thought of as obstacles for the formation of micro-crystallites) should result in 
a less pronounced increase of the average time for reformation of strands, r^,, with the 
growth of macro-strains. This conclusion is confirmed by the data listed in Table 2 which 
reveal that the parameter ri in Eq. (59) (the "rate of growth" for the characteristic time 
of relaxation with the measure of deformation, Ji — 3) for the unfilled rubber exceeds that 
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for the filled rubbers by several times. 

An increase in the average relaxation time, r*, with stretching of a specimen is correctly 
predicted by Eq. (59) for the unfilled rubber and for the filled rubber Rl. For the other 
two grades of particle-reinforced rubber (R2 and R3), T^, falls down rapidly in the interval 
between A = 1.4 and A = 1.8 and proceeds to increase for larger values of the elongation 
ratio. This decrease in (or the corresponding increase in the rate of breakage F^) 
may be explained if we assume that links between filler particles and long chains are 
strong enough in order to prevent thermally induced detachment of strands from the 
particles at relatively small strains. This implies that in the region of relatively small 
elongation ratios (A < 1.4) particles of filler may be treated as permanent cross-links 
between chains. When the macro-deformation of a specimen reaches a threshold value for 
the load-bearing capacity of these links, chains begin to slip from the particles, which is 
reflected by the model as a noticeable increase in the rate of breakage F*. The parameter 
Ti (that describes the effect of stretching on the characteristic time of relaxation) in the 
region of low elongation ratios (0 < A < 1.4) substantially exceeds that in the domain 
of large stretches (1.8 < A < 2.0). This may be explained by the interaction of two 
physical phenomena in the region of large deformations: (i) the growth of A results in 
slowing down of the reformation rate (which is in common for regions of low and high 
strains), and (ii) an increase in the concentration of temporary cross-links (when the 
stretch intensity exceeds the load-bearing capacity of links between macromolecules and 
filler particles) accelerates the reformation process (which is a specific feature for the 
domain of relatively large strains). 

It is worth noting that the parameter ti for the filled rubber R2 (samples cut out in 
the direction orthogonal to the direction of milling) noticeably exceeds that for the filled 
rubber Rl (specimens cut in the direction of milling). This may be explained by the 
mechanically induced crystallization of specimens before their vulcanization. Chains with 
the end-to-end vectors directed along the direction of loading in the relaxation tests are 
partially crystallized during milling, which results in an essential decrease in their rate 
of breakage and reformation. Because these chains bear the main load under tension for 
the specimens Rl, the influence of macro-deformation on the breakage process becomes 
relatively weak (small values of Ti). On the contrary, partial crystallization of chains with 
the end-to-end vectors directed perpendicular to the direction of stretching for specimens 
R2 does not affect their rates of breakage (because micro-strains in these chains are 
relatively small), whereas the chains directed along the stretching direction are practically 
not crystallized, and their response is similar to that for the unfilled rubber (large values 

ofTi). 

Figures 3 and 7 demonstrate that A increases with the longitudinal strain, and this 
dependence can be correctly described by the phenomenological equation (59). According 
to Table 2, the value of A in the stress-free state, Aq, equals approximately 0.2 for all 
compounds under consideration. Bearing in mind that A is the ratio of the stress in 
regions with low concentration of junctions to the total stress in a specimen, see Eq. (58), 
and presuming (as a first approximation) that the stress is linearly proportional to the 
mechanical energy, we conclude that the mechanical energy for RLCJs equals about 0.2 of 
the total strain energy for rubber. Finally, assuming that the energy stored in a domain 
is proportional to its volume, we arrive at the rough estimate for the volume occupied by 
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regions where relaxation of stresses occurs: about 20 % of the total volume of a specimen 
(Table 2 shows that this value weakly depends on the content of filler). 

Equation (59) implies that an increase in the parameter A with the elongation ratio 
A means that the stress, ctl, in RLCJs grows more rapidly then the stresses in the bulk 
material, ctb, and in the regions with high concentration of junctions, cth- This conclusion 
is fulfilled both for the unfilled rubber and for the filled rubbers Rl and R2 (which were 
subjected to milling during the preparation process). On the contrary. Figure 10 demon- 
strates that at relatively small elongation ratios, A < 1.6, the parameter A decreases, 
whereas at large stretches, A > 1.7, it linearly increases with the first principal invariant, 
/i, of the Cauchy deformation tensor, C°. This non-monotonicity in the dependence of 
A on strains for specimens R3 may be explained within the concept of occluded rubber 
38| . According to that theory, aggregation of isolated clusters of particles into a network 



results in the separation of some regions of the host material (occluded domains) from 
the bulk elastomer. When the secondary network of filler is relatively rigid, it does not 
transmit stresses to the occluded domains. As a consequence, at small strains, occluded 
rubber is not deformed, and its mechanical energy vanishes. With the growth of strains, 
the network of filler breaks, and some regions of the occluded rubber are released from the 
constrains. This results in an increase in the total strain energy and stresses in the bulk 
material. With reference to Eq. (58), this conclusion implies that in the region of rela- 
tively small elongation ratios (A < 1.6) for rubber R3, the sum ctb + cth grows faster then 
the stress, ctl, in the domains with low concentration of junctions. When the elongation 
ratio A reaches its threshold value (A ~ 1.7), all regions of occluded rubber are released 
from restrains, and the growth of strains implies a more rapid increase in the stress ctl 
then in the stresses in the host elastomer and in the regions with high concentration of 
junctions. The latter is observed as the growth in the parameter A. The difference in the 
behavior of the function A{Ii) for the unfilled rubber and for the CB filled rubbers Rl 
and R2, on the one hand, and for the rubber R3, on the other hand, may be explained 
by the absence of occluded regions for the first three compounds (the regions of occluded 
rubber are supposed to be destroyed during milling) and the presence of these regions in 
the samples R3 with no preliminary mechanical treatment. 

It follows from Figures 3 and 7 that for any elongation ratio A > 0, the quantity A 
is maximal for the unfilled rubber, and the value of A for specimens cut in the direction 
orthogonal to the direction of milling (R2) is higher then that for specimens cut in the 
direction of milling (Rl). The fact that the ratio A for the unfilled rubber exceeds that 
for the filled rubbers (Rl and R2) follows from Eq. (58), provided that the stress in 
the regions with high concentration of junctions exceeds that for the bulk material. The 
observation that the parameter A for specimens R2 exceeds that for the specimens Rl 
confirms our hypothesis about the effect of mechanically induced crystallization on the 
response of filled rubbers, because ub increases when macromolecules whose end-to-end 
vectors are directed along the axis of extension are partially crystallized. 
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8 Concluding remarks 



Constitutive equations have been derived for the isothermal viscoelastic behavior of un- 
filled and filled elastomers at finite strains. The model is based on the concept of tem- 
porary networks that treats a rubbery polymer as a network of long chains connected 
to junctions (chemical and physical cross-links and entanglements). A filled rubber is 
thought of as a strongly non-homogeneous medium (with the characteristic length of in- 
homogeneity of about 100 nm) which consists of an elastomeric matrix with a normal 
concentration of cross-links, wherein regions with low and high concentrations of cross- 
links are distributed. In the domains with low concentration of cross-links, the average 
number of junctions between polymeric chains (per unit mass) is less then that in the 
bulk material because of local inhomogeneities in the distribution of a cross-linker during 
the mixing procedure. In the domains with high concentration of cross-links, the number 
of junctions between long chains exceeds that for the host elastomer, because particles of 
filler and their clusters play the role of extra cross-links. 

In the bulk material and in the regions with high concentration of cross-links, inter- 
molecular forces prevent slippage of chains from junctions (which are treated as permanent 
cross-links). In the domains with low concentration of cross-links, active strands (whose 
ends are linked to separate junctions) are assumed to slip from the temporary junctions as 
they are thermally agitated. Breakage of an active strand (its transition into a dangling 
strand) occurs at a random instant when one of its ends detaches from a junction. A 
dangling strand is reformed (it merges with the network) at a random time when its free 
end captures some junction in its vicinity. The stress in a dangling strand is assumed to 
totally relax before this strand is attached to the network, which implies that the natural 
(stress-free) state of any chain coincides with the deformed state of the network at the 
instant of their merging. Unlike conventional theories of transient networks, we postulate 
that different strands in RLCJs have different energies for breakage and describe these 
energies by the quasi-Gaussian distribution function (12) with two adjustable parameters, 
Q and S. 

Constitutive equations for a filled rubber are developed by using the laws of thermody- 
namics. These equations are applied to determine stresses at uniaxial extension of a bar. 
The time-dependent response of a specimen in the standard relaxation test is character- 
ized by four adjustable parameters which are found by fitting observations. To determine 
these constants, a series of tests have been carried out on specimens of unfilled and CB 
filled natural rubber. Fair agreement is demonstrated between results of numerical simu- 
lation and observations in relaxation tests (with the duration of to = 1 hour) at various 
elongation ratios. A, in the region of moderate finite deformations, A e (1, 2). 

The following conclusions are drawn: 

1. The unfilled and particle reinforced rubbers demonstrate a pronounced viscoelastic 
response. The relative decrease in stresses, C(^o) , for the filled rubbers exceeds that 
for the unfilled rubber by twice. This difference may be ascribed to an increase 
in the inhomogeneity of spatial distribution of junctions in an elastomeric matrix 
caused by the presence of filler particles. 

2. The characteristic time of relaxation, r*, for the unfilled and filled rubbers increases 
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with the elongation ratio, A. This observation contradicts the strain-time super- 
position principle conventionally employed for the analysis of glassy polymers. An 
increase in r* with A may be explained by the mechanically induced crystallization 
of natural rubber in the regions with low concentration of junctions. The rate of 
growth for the average relaxation time, r*, with strains is maximal for the unfilled 
rubber and it decreases with an increase in the content of filler. 

3. The average energy for breakage of strands, Q, decreases with an increase in the 
filler content, whereas the width of the Gaussian distribution of energies, ^ = S/f2, 
is practically independent of the concentration of carbon black. The distribution 
function for the energies of breakage, p{uj), does not depend on the elongation ratio, 
A. 

4. The ratio of stresses in the regions with low concentration of junctions to the total 
stress in a specimen. A, increases with A. This means that at uniaxial tension, 
stresses in the domains with low concentration of cross-links grow more rapidly then 
in the elastomeric matrix and in the regions with high concentration of junctions. 
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Table 1: Chemical composition of rubbers (phr by weight) used in the experimental 

study 



Formulation 


EDS-16 


EDS-19 


Natural rubber 


100 


100 


Zinc oxide 


5 


5 


Stearic acid 


2 


2 


Carbon black (N330) 


45 


1 


Process oil 


4.5 





Antidegradant(HPPD) 


3 


3 


Antiozonant wax 


2 


2 


Accelerator (CBS) 


0.6 


0.6 


Sulphur 


2.5 


2.5 



Table 2: Adjustable parameters of the model 



Rubber 


W 


S 




^0 


Ai 


To 




Unfilled 


11.70 


6.10 


0.52 


0.20 


0.054 


0.00 


9.20 


Rl 


5.30 


2.80 


0.53 


0.17 


0.024 


1.01 


0.12 


R2 


5.13 


3.18 


0.62 


0.18 


0.026 


0.68/1.30 


3.02/0.30 


R3 


5.55 


3.36 


0.61 


0.15/0.07 


-0.07/0.02 


0.54/0.22 


2.78/0.89 
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List of figures 



Figure 1: The dimensionless ratio R versus time t (s) for unfilled rubber in relaxation 
tests with the elongation ratios A at room temperature. Symbols: experimental data. 
Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; curve 3: 
A = 1.8 

Figure 2: The characteristic time of relaxation r^, (s) versus the first invariant of the 
Cauchy deformation tensor Ii for unfilled rubber in relaxation tests at room temperature. 
Circles: treatment of observations. Solid line: approximation of experimental data by Eq. 
(59) with To = and n = 9.2012 

Figure 3: The dimensionless parameter A versus the first invariant of the Cauchy de- 
formation tensor Ji for unfilled rubber in relaxation tests at room temperature. Circles: 
treatment of observations. Solid line: approximation of experimental data by Eq. (59) 
with Ao = 0.2028 and Ai = 0.0537 

Figure 4: The dimensionless ratio R versus time t (s) for CB filled rubber (Rl) in 
relaxation tests with the elongation ratios A at room temperature. Circles: experimental 
data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; 
curve 3: A = 1.8; curve 4: A = 2.0 

Figure 5: The dimensionless ratio R versus time t (s) for CB filled rubber (R2) in 
relaxation tests with the elongation ratios A at room temperature. Circles: experimental 
data. Solid lines: results of numerical simulation. Curve 1: A = 1.2; curve 2: A = 1.4; 
curve 3: A = 1.8; curve 4: A = 2.0 

Figure 6: The characteristic time of relaxation r* (s) versus the first invariant Ji of the 
Cauchy deformation tensor for CB filled rubbers in relaxation tests at room temperature. 
Symbols: treatment of observations. Unfilled circles: Rl; filled circles: R2. Solid lines: 
approximation of experimental data by Eq. (59). Curve 1: tq = 1.0114, ri = 0.1217; 
curve 2a: tq = 0.6777, n = 3.0218; curve 2b: tq = 1.2954, n = 0.2957 
Figure 7: The dimensionless parameter A versus the first invariant of the Cauchy de- 
formation tensor Ji for two CB filled rubbers in relaxation tests at room temperature. 
Symbols: treatment of observations. Unfilled circles: Rl; filled circles: R2. Solid lines: 
approximation of experimental data by Eq. (59). Curve 1: Aq = 0.1677, Ai = 0.0237; 
curve 2: Aq = 0.1762, Ai = 0.0263 

Figure 8: The dimensionless ratio R versus time t (s) for CB filled rubber in relaxation 
tests with the elongation ratios A at room temperature. Circles: experimental data |2T|. 
Solid lines: results of numerical simulation. Curve 1: A = 1.19; curve 2: A = 1.38; curve 
3: A = 1.58; curve 4: A = 1.76; curve 5: A = 1.92 

Figure 9: The characteristic time of relaxation r* (s) versus the first invariant of the 
Cauchy deformation tensor Ji in relaxation tests at room temperature. Circles: treatment 
of observations |]21|. Solid lines: approximation of the experimental data by Eq. (59). 



Curve 1: tq = 0.5373, n = 2.7795; curve 2: tq = 0.2204, n = 0.8877 
Figure 10: The dimensionless parameter A versus the first invariant of the Cauchy 
deformation tensor Ji for relaxation tests at room temperature. Circles: treatment of 
observations [PT|. Solid lines: approximation of the experimental data by Eq. (59). 
Curve 1: Aq = 0.1467, Ai = -0.0702; curve 2: Aq = 0.0697, Ai = 0.0206 
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